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Abstract. We compute the three-point cross-correlation function of the primordial curva- 
ture perturbation generated during inflation with two powers of a vector field in a model 
where conformal invariance is broken by a direct coupling of the vector field with the in- 
flaton. If the vector field is identified with the electromagnetic field, this correlation would 
be a non-Gaussian signature of primordial magnetic fields generated during inflation. We 
find that the signal is maximized for the flattened configuration where the wave number of 
the curvature perturbation is twice that of the vector field and in this limit, the magnetic 
non- linear parameter becomes as large as \bNL\ ~ 0(10^). In the squeezed limit where the 
wave number of the curvature perturbation vanishes, our results agree with the magnetic 
consistency relation derived in arXiv:1207.4187. 
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1 Introduction 

The observations of the acoustic peaks in the Cosmic Microwave Background (CMB) and the 
underlying nearly flat spectrum of primordial curvature perturbations are so far the main 
observational veriflcation of the inflationary paradigm and has ruled out some of the leading 
competing scenarios such as cosmic string structure formation [1]. Despite the success of 
the inflationary paradigm any flrm knowledge about the underlying microscopic model is 
still lacking. In fact, even the energy scale of inflation is still unknown within ten orders of 
magnitude, but a future positive detection of primordial gravitational waves would be able 
to pin it down. Similarly a detection of primordial non-Gaussianity or isocurvature modes 
would give us more detailed information about the underlying model of inflation. These are 
some of the new observables, that we might be fortunate to probe with Planck and other 
present and future experiments, but it is strongly model dependent if these observables will 
turn out to be within experimental reach [2]. It is therefore important to keep looking for 
new observables, which could be significant in other non-minimal versions of infiation. 

On the other hand astrophysical observations indicate the existence of coherent magnetic 
fields on different cosmological scales. The strength of magnetic fields in galaxies and galaxy 
clusters has been determined to be of the order of a micro Gauss [3], while a lower bound 
of 0.1 femto Gauss on magnetic fields in the intergalactic medium has also been claimed [4]. 
While it may be hard to see how astrophysical processes can generate coherent magnetic fields 
on such large scales, inflation by its very nature creates coherent fluctuations on very large 
scales. It is therefore natural to speculate that large scale primordial magnetic fields could 
be generated during inflation by a coupling between the inflaton field and the vector field, 
which breaks the conformal invariance of the vector field during inflation [5]. As we discuss 
below, such models of primordial magnetogenesis suffer from the strong coupling problem [6] . 
However, if gauge flelds (which may or may not be identified with the electromagnetic field) 
are produced this way, they will be correlated with the primordial curvature perturbation 
generated during inflation. 

Motivated by these observations, the focus of the present paper is therefore to compute 
the three-point cross-correlation function between the comoving curvature perturbation and 
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two powers of the gauge field (to be identified with the electromagnetic field in models of 
magnetogenesis). We have recently derived a magnetic consistency relation [7], which can be 
used to check the results in the so-called squeezed limit, and we find very good agreement 
between the consistency relation and the results presented here. 

Related calculations, which have partly served, as motivation for the present work, 
have previously been carried out in [8] in the more restricted setting of a test scalar field 
e.g. curvaton [9-11] coupled to electromagnetism. Subsequently, the cross-correlation of 
the comoving curvature perturbation with the magnetic fields in a specific model was also 
considered [12]. Our calculations can be regarded as a generalization and check of these 
results, although there are some disagreements, as we will discuss later. 

The time-dependent coupling of the U{1) gauge field to the background (or the scalar 
field) can be parametrized by a coupling of the form A((/))F^,yF^'^, where F^^, = d^j^Ay — d^A^ 
is the field strength associated to the U{1) gauge field, and the time dependence of the 
coupling is parametrized by the dependence on some slowly rolling background field (f), which 
we will think of as being the inflaton for simplicity. If the vector field is identified with the 
electromagnetic field this is a well studied model of primordial magnetogenesis [5, 13, 14]. It 
has been proposed that such type of models could yield the required magnetic seed fields for 
the galactic dynamo [15-23], and perhaps lead to signatures in the CMB [24-38]. 

We would like to emphasize that there are some theoretical obstacles for a model of this 
type to succeed as a working scenario of primordial magnetogenesis, when the gauge field 
is associated with the electromagnetic field. In order to produce significant magnetic fields 
without running into problems with large back reaction on large or small scales, the spectrum 
of the magnetic fields has to be nearly scale invariant. On the other hand, a scale invariant 
spectrum of magnetic fields can only be obtained for n = 2 and n = —3, when the coupling 
A((/>) ~ a'^{t) with a{t) being the scale factor, varies by many orders of magnitude during 
inflation. This implies that for n = 2, if \{4') at the end of inflation is unity, i.e.. A/ = 1, as 
required to obtain standard electromagnetism, the effective gauge coupling squared, scaling 
inversely with A, must have been very large at the beginning of inflation and one would 
have started out in a regime of strong coupling [6]. On the other hand, when n = — 3, 
which is in the weakly coupled regime, it has been shown that the energy density in the 
electromagnetic field scales as £em ~ iJf a"^^"'"'"^^ . Therefore, for backreaction to remain 
small, i.e., £em ^ H^-, for a typical value of Hj ~ 10"^, we can maximally obtain only about 
12 e-folds of infiation before backreaction becomes important, which is much less than the 
required number of e-folds. 

As discussed in [39] it is very difficult to see how the strong coupling issue can be avoided 
retaining gauge invariance, unless one gives up on the relation to primordial magnetogenesis 
and on the gauge field being the actual electromagnetic field, and instead thinks of it as some 
novel hidden sector isocurvature gauge field. Another possibility, which has been proposed, 
is that the gauge invariance is broken in the UV in the effective four dimensional Lagrangian, 
but manifest in five dimensions [12]. If gauge invariance is broken, the effective gauge coupling 
squared may not scale inversely with It has also been proposed that the gauge field 

could instead be identified with a vector curvaton, in which case the constraint on the gauge 
coupling at the end of infiation doesn't apply. In this case we would have a new type of 
non-Gaussianity in a scenario of mixed infiation and vector curvaton perturbations [7]. Here 
we will be agnostic about the detailed role of the gauge field, and leave further exploration of 
possible applications of the model for future work. However, for the remaining of this paper, 
we will for definiteness refer to the gauge field as if it is the electromagnetic field, which will 
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also make the comparison with the previous work in [8, 12] more direct. 

It is not immediately obvious what is the most convenient way to parametrize the 
results. If we define the cross-correlation bispectrum of the curvature perturbation with the 
magnetic fields as 

(C(ki)B(k2) • B(k3)) = (27r)35(3)(ki + k2 + k3)i?CBB(ki,k2,k3) , (1.1) 

then as we have previously proposed, it is convenient to define the magnetic non-linearity 
parameter^ ^tvl, in terms of the cross-correlation function of the curvature perturbation with 
the magnetic fields 

S^BB(ki,k2,k3) = bNLPc{h)PB{k2) , (1.2) 

where and Pb are the power spectra of the comoving curvature perturbation and the 
magnetic fields, defined respectively as 

(C(k)C(k')> ^ (2^)35(3)(k + k')Pc(^), (1-3) 
(B(k)-B(k')> = (27r)35(3)(k + k')PB(fe). (1.4) 

In the case where b^i is momentum independent, it takes a "local" form which can be 
derived from the relation 

B = b(«) + h^f^C^'^^Bi^^ (1.5) 

with B^*^) and C^'^^ being the Gaussian fields. There is an interesting limit where the magnetic 
non-linearity parameter takes the local form, which makes the comparison with the above 
estimate particularly simple. We will show that in the squeezed limit, where the momentum 
of the curvature perturbation vanishes, i.e., ki <C ^2,^:3 = k, we, in fact, recover 

(C(ki)B(k2)-B(k3)) = 6^^(2^)35(3) (k,+k2 + k3)Pc(^l)^i?(^) , (1-6) 

with 6^™' = ub — 4: where is the spectral index of the magnetic field power spectrum, 
in agreement with the magnetic consistency relation, which was derived in [7] using simpler 
semi-classical methods^. In the case of a scale invariant spectrum of magnetic fields, ub = 0, 
we have ^Jv*^' = — 4. 

Another interesting limit which maximizes the three-point cross-correlation function is 
the fiattened shape where ki/2 = ^2 = ^3. In this limit it turns out that the signal is 
enhanced by a logarithmic factor in agreement with [8]. On the largest scales the logarithm 
will give an enhancement by a factor 60. Thus, for a flat magnetic field power spectrum, 
the non-linearity parameter in the fiattened limit becomes I^tvlI ~ 0(10^) depending on the 
scale. 

The paper is organized as follows. In the next section we will review the model and 
the mechanism for the production of large scale gauge fields during inflation, which has been 
proposed as a mechanism of primordial magnetogenesis. In section 3, we derive the interaction 
Hamiltonian for the primordial curvature perturbation and the gauge fleld, while in section 
4 we calculate the three-point cross-correlation function of the curvature perturbation with 
two powers of the magnetic field. We show that ^at^ is maximal in the flattened shape and 

^Note, a different dimensionless quantity Bqbb (ki, k2, / PQ{k'i)PB{k2)PB{k-3)) was previously intro- 
duced in [8, 12], which is different from our bjvL 

^This approach is a non-trivial generalization of related semiclassical methods used in [40-46]. In the 
appendix of [43], some of these approaches are reviewed. 
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of order O(IO^). We then derive the squeezed hmit of our result, and compare it with the 
magnetic consistency relation. Finally, in section 5, we conclude with a discussion of our 
results. In the appendix we have listed some useful integrals. 

Throughout this paper, we work in natural units with h = c = 1, and the Planck mass 
Mp = l/SvrG equal to unity. Our metric convention is (—,+,+,+). 

2 Primordial magnetogenesis 

In order to break the conformal invariance of the Maxwell equations, a necessary condition 
for amplification of magnetic fields during inflation, we consider the coupling of the electro- 
magnetic field with a scalar field [5] 

^=\J d^x^ (r - d^<Pd^^ - '^y{<P) - ^A(</.)F^,F'^^^ . (2.1) 

In general, the scalar field doesn't need to be the inflaton, but could be any light field (or 
a set of light fields) in slow-roll during infiation. However, for simplicity we will assume by 
default that (p is the infiaton, and comment in the appropriate places how the result would 
change if (f) is some other light field like a curvaton or an entropy perturbation mode. If is 
the inflaton field, it has to satisfy the usual slow-roll conditions 

1 fV'V V" 

We find it convenient to work in the Coulomb gauge with = and diA^ = 0, where 
the quadratic action for the electromagnetic vector field, Aj, becomes 

Sem = -\j d'x^A(0)F^,F^'^ = ^ y SxdTX{(l)) (^Af - ^{drAj - djAif^ . (2.3) 

The conformal time, r, and the scale factor, a(T), are defined by the choice of the background 
metric ds^ = a^(T)(-(ir^ + dx^). 

We are interested in calculating correlation functions in Fourier space as in (1.2), where 
the Fourier transformed gauge field is 

Ai(T,k) = j d''xAi{T,i^y^-'' (2.4) 

In order to quantize the gauge field, we define the usual mode expansion 

Ai{T, k) = [e1{k)Ak{T)bl^''-^ + h.c] (2.5) 

(T=± 

and impose the standard commutation relations 

[6£,6-;V(2vr)=^5(3)(k-k')5..' . (2.6) 
The polarization vectors will have to satisfy 

k • e%k) = , e%k) ■ e^'*(A:) = 5^^, , e1{k)e^^*{k) = <5,, - kik./k' . (2.7) 

CT=± 
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For the two point correlation function, we then have 

{A{T,k)AjiT,k')) = (27r)3d(3)(k + k') (^s,, - ^) |A,(r)|2 . (2.8) 

We can now compute the magnetic field in terms of the gauge field, by using the relation 

Bi{T,x) = ^eijkdjAk{T,x) . (2.9) 

Upon using eij^e*™" = ^J^^k ~ ^^^k^^ that the two point correlation function of the 

magnetic fields is given by 

= 2^(27r)35(3)(k + k')|ylfc(r)|2 (2.10) 

If we define the magnetic power spectrum as 

{B,{T,k)B\T,k')) = {27rfS^^\k + k')PBik) , (2.11) 
then the magnetic power spectrum is given by 



PB{k) = 2^\Ak{T)\\ (2.12) 



It is convenient to parametrize the time dependence of the coupling as A((/>(r)) = 
Xi{t/tj)~'^^. Defining the pump field S'^{r]) = X{(j){r])) and the canonically normalized vector 
field Vi = S{T)Ai, the quadratic action for the vector field takes the canonical form for the 
kinetic term 

Sem = 2 ^ drd^x 



S 

and the equation of motion for the mode function, = S{T)Ak, is 



(2.13) 



v'l+{k'-^]vk = 0. (2.14) 



The solution, normalized to the Bunch-Davies vacuum, is 



Vk{r) = :^e^-('^^>^^V^H'^lJ-kT) , (2.15) 



which leads to [8] 



Psik) = j-^l^ i-krfH^^li-kr)Hfli-kr) (2.16) 

Assuming a resolution to the strong coupling problem as in [12], it can be estimated that for 
n = 2, where backreaction remains small and with H ~ 10^^ GeV, a magnetic field strength 
of ~ 1 nano Gauss can be achieved on Mpc scales [6, 8]. 
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3 The interaction Hamiltonian 



To find the interaction Hamiltonian of tlie metric fluctuations with the electromagnetic fields, 
it is convenient to write the metric in the ADM form as 



N'^dt^ + hij{dx' + N'dt){dx^ + N^dt) , 



(3.1) 



where the dynamic degrees of freedom are contained in hij, and the lapse, N, and the shift, 
N^, are determined by the constraint equations, which can be derived from the Lagrangian 
obtained by inserting the ADM decomposition of the metric into the background action 



iVi?(3) - 2NA + N^'^iElE] - {Elf) 



2 

+ N~'^{^-N'di 



?-Na~^\e 



71«J 



(3.2) 



where Eij = ^{hij —\7 iNj —\7 jNi) is the rescaled extrinsic curvature and R^^^ is the curvature 



scalar of the three-metric hij. 



The dynamical degrees of freedom can be parametrized by a scalar curvature perturba- 
tion, and traceless and transverse tensor perturbation, 7jj, by writing 



hi 



In this way, one obtains for the lapse and the shift to linear order 



(3.3) 



N 

Ni = a 



1 

H 



Similarly quadratic action for the curvature perturbation becomes 



S 



d^xdTa\T)e[C'^ -{dCf] . 



Expanding in terms of mode functions 

C(r,x) = 
the mode equation becomes 



d^k r 
(2^ 



Q'l + 2aH{l + 6 + e)C', + = 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



with e = —H/H^ and 6 = H/2HH. The solution normalized to the Bunch-Davis vacuum to 
leading order in slow-roll is given by 



Ck{r) 



1 H 



(1 -|- ikr) e' 



-ikr 



(3.8) 



Under linear variation of the Lagrangian in the metric, the interaction Hamiltonian is 
given by 

HcAA = -^jd^xa^ T^'''^g^, . (3.9) 
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Upon using (3.4) in (3.9), the interactions of the metric fluctuations with gauge fields at third 
order are given by [47] 

HcAA = -\j d'xa' (Icroo - a. (-Ic + ahd-'C^ T"^ - a\T^^ . (3.10) 

Inserting directly the energy-momentum tensor of the gauge field 

= AF^pF/ - ^g^.XFp^F"'' , (3.11) 

one finds an expression that agrees with equation (46) of [12] and is of order \(AA. However, 
the true interaction is actually of lower order since the interaction used in [12] is proportional 
to a total time derivative to leading order. Using that 



V^T'^'^ = -iVX)Fp^FP^ , (3.12) 



which in the comoving Coulumb gauge yields 

^f^'^^'' = -^>^[^^^—^i9iAj-djA,f^ , (3.13) 

where the non-vanishing right-hand-side is due to the sourcing between the salad and vector 
modes. The total energy-momentum tensor of the gauge field and the inflaton is of course 
conserved. 

Writing explicitly the covariant derivative in terms of the background metric, one finds 
^dt (a^rOO) + aaT" + 9,^° + (^i^. - ^(diA, - d.Aif^ = , (3.14) 

so clearly, by inserting (3.14) into (3.10), we obtain after some partial integrations to leading 
order in slow-roll 

HcAA = -\j d'x {aX^C {a,A, - ^{^^Aj - d,A,f^ - dt («'|:CT°°)) . (3.15) 

In the above expression, we have dropped a total spatial derivative term. The total time 
derivative can be removed by a field redefinition, which will however only affect the results 
at higher order in perturbation theory. Using 6(f) = \/2eC ^iid -v/2e = —(j)/H, we have 

and with this substitution, the action agrees to leading order in slow-roll with the action in 
the uniform curvature gauge given in equation (2.23) of [48] and with equation (30) of [8] in 
the special case discussed there. 

Note that this simple form of the action enables us to generalize our results, in a 
straightforward manner, to the case where the scalar field 4> is not the inflaton but either 
an isocurvature field or a curvaton, since the action will take the same simple form in these 
cases. 
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4 The correlation of curvature perturbations with the magnetic fields 



In order to compute the higher-order correlation function during inflation, we adopt a very 
useful and powerful tool of the in-in formalism [40]. In this formalism, the expectation value 
of an operator O at time r/ is given by 



(0| OiTi) \n) = (0| T (e'f-'^'^^^A OiTi)T fe-*/-- 



drH, 



Olo) 



(4.1) 

where |0) is the vacuum of the interacting theory, |0) is the vacuum of the free theory, T 
and T are time ordering and anti-ordering operators, respectively, and Hj is the interaction 
Hamiltonian for time r. 

From the interaction Hamiltonian in (3.15), we obtain the cross-correlation of the cur- 
vature perturbation with the electromagnetic field 

{C{Tl,ki)MriM)A,{Tj,k3)) = (27r)35(3)(ki+k2 + k3) 



Si 



^2 



Oil 79— I ksi Ojk p5 ) K2,fc-i^2 



^2 



with the integrals 



23= 



CkAri)AkAri)Ak.Ari) / drrA'(r)a,(T)9. (^^.^(r)) 9. (^^^(r)) 



and 



23= 



CkAri)Ak,{Ti)Ak,{Ti) / dTTX'iT)aM)Al{T)Al^{T) 



(4.2) 



(4.3) 



(4.4) 



where the mode function CkiT) is given in (3.8) and the mode function ^^(t) is obtained 
from (2.15) which is given by 



Akir) 



1 



,j7r(l+n)/2 



A/ 2 



(1) 



-kr) 



(4.5) 



Inserting the expressions for the mode functions, and using the property of the Hankel func- 
tions 

dJx-Hi'''Hx))=x-H^^l'^{x) , (4.6) 



we can rewrite the two integrals as 

Xi = (-2n) |cIf(r,)|2|4?(T,)||4°)(r,)| ^2^3 



2n) lC(r,)|Vi?(r/)|l4?(-/)l^l^^ 



(0), 



(4.7) 
(4.8) 



where Ck^\Tj) and A^^\tj) are the asymptotic super horizon values of the mode functions 
and are given by 



lCf(r.)| 



1 H 

MV2k^ ' 

1 r(n+l/2)2„_i 1 



AT r(3/2) 



2k 



-kriY 



(4.9) 
(4.10) 
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and the integrals are 



-2n-l 



2 r2(n+ 1/2) 



(1 + .A;ir/)e"^^-^-^ifi'_),/2(-A;2r,)i?(i;i/2(-^3r/) 



^ii/2(-^3r) 



(4.11) 



.(_fe,^,)"+l/2(_fc3^^)n+l/2 



2 r2(n + 1/2) 

a + ^A:ir,)e-^'=i-^ili^_),/2(-A;2r/)i?i!;i/2(-^3r/) 

' drr(l - ifeir)e^^i^/?^'^^/2(-^2T)F 



^(2) . 
n+1/2^ 



/csr) 



(4.12) 



Using dlna = Hdt, one has — 2n 
correlation function becomes 



-(aA//aina)/A7 



-Xi/{HXj), and in this way the 



{C{TjM)A{TjM)MTj,ks)) = -l^(2vr)35(3)(k, + k2 + k3)|Cif(T/)P|4?(T,)||4°)(T,) 



k2,ik2,i 

^2 



5ii 



j k3i [djk 



^ ) k2,A'^ 

K/Q 



(4.13) 



This expression represents our most general result, where the gauge field indices have not 
yet been contracted. The cross-correlation of the curvature perturbation with the magnetic 
fields is now given by 

{({tiMMtiM) ■ B(T/,k3)) = {6,jk2 ■ ks - k2,,k3j) (C(r/ , ki)^, (tj, ks) A, (t/ , ks)) 



(4.14) 



which leads to the final result 



(C(r,,ki)B(T,,k2)-B(r,,k3)) = l^(2vr)35(3)(ki + k^ + k3)|ci;^(r,)|2|4?(T,)||4°^(r, 



|2| /l(O), 



1(0), 



HXi 



k2 • ks + 



(k2 • k3)3 

«,2«'3 



k2ksii^^ + 2(k2 • k3)2x(2) 



(4.15) 



Note that the polarization factor multiplying the integrals is slightly different than in equation 
(41) of [8], where it appears that part of the polarization tensor product was missed. 

The solutions of the integrals for different values of n are listed in the appendix. How- 
ever, for the most interesting case of n = 2, we find 



1 



{k2k3)y^kf 



(4.16) 



[-kf - 2kl{k2 + k^) - 2ki{kl + A;2A:3 + kl) - (^2 + k^){kl + ^2^:3 + kl)] 
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and 



X [(A;i + k2f{-3kl - 3/c?/c2 - k^) + {ki + k2){-9kl - 6klk2 - 2k^)k3 
+{-9kf - 6klk2 - 2kikl - 2kl)kl 

-2{2kl + kik2 + kl)kl - 2{ki + k2)k^ - 4 + 3kfkf{j + ln(-A;tT7))] 
where we have defined kt = ki +k2 + k^ and 7 is the Euler gamma constant. 
4.1 The flattened shape 

It is interesting to note that the ln(— /c^r/) term only appears in (4.17), and can therefore 
not cancel out in general. This term will be most important when ki is maximized in the 
flattened shape with ki = 2k2 = 2ks as already observed in [8]. For length scales relevant for 
CMB, the logarithm will give an enhancement by a factor 60, but on smaller scales it can be 
even larger. Since the logarithm completely dominates the integral in the flattened limit, it 
is easy to estimate the size of h^i in this limit. When the logarithmic term dominates, we 
have 

The contribution of this term to the non-Gaussian cross-correlation function is 

(C(T,,kOB(r,,k2)-B(r,,k3)):^6-i^ ^^^^^'^^^^ ln(-fc,r,)|cif (rj)| Vg(T7)|l4?(^/)l • 

(4.19) 

For n = 2, we have {Xj/HXj) = —2n = —4, and using ki = 2k2 = 2hi in the flattened limit, 
we have 

(C(r/, ki)B(r7,k2) • B(rj,k3)) ~ -m\n{-ktTi)P^{ki)PB{k2) ■ (4.20) 

For the reasonable value of the logarithm, corresponding to the largest observable scale today, 
\a{—ktTi) ~ 60, we then obtain in the flattened limit 



,flat 



5760 . (4.21) 



This is a quite signiflcant contribution to the non-Gaussianity as compared to the a priori 
expected level of order unity. 

4.2 The squeezed limit 

Now let us consider the special limit where the wavelength of the curvature perturbation is 
much longer than the wavelength of the magnetic fields. In this limit we have ki — t- and 
k3 — )• — k2 = — k. Using the asymptotic behavior for the real and imaginary parts of the 
Hankel function ^[Hi^\x)] and "^siH^Pix)] cc X ^ for X —7- 0, it is possible to verify, 

that in the squeezed limit, the integrals reduce to 

X« = j' dTJn^y2{-kTi)Y^-i/2{-kTi) (4.22) 

and 

ii?'=i!di- (4.23) 
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In the squeezed limit, the cross-correlation therefore reduces to 



X 



(4-24) 



For integer values of n, it can be proven that the integral in (4.22) is 

2<i) = -(n- 1/2)^2 (4.25) 

which (in the ki — )• limit) gives 

(C(r/,ki)^,(T7,k2)A,(T/,k3)) = --^^(27r)V3)(ki + k2 + k3) 



x(^.-i^)icr(-^)ivr(r,)p 



iCiri, ki)C(T7, -ki)) {Mn, k2)MTl, k3))(4.26) 



We have shown numerically that (4.25) also holds for real non-integer values of n. 

For the cross-correlation of the curvature perturbations with the magnetic fields, we 
then obtain for n > 

(C(t/, ki)B(r7, k2) • B(r/, k3)) =--|^(27r)35(3)(ki + k2 + k3)Pc(A;i)^'B(A;2) • (4.27) 

This agrees with the squeezed limit result in equation (64) of [8], when using AC = —d^\H5(j) 
and inserting the specific form of the coupling = exp(2(f>/M) used there. This agreement 
is however a coincidence because the difference in the polarization sums noted after (4.15) 
vanishes in the squeezed limit. In fact, if we had not taken the trace of BiBj in the correlation 
function, the results would no longer agree, even in the squeezed limit. We also note that both 
of these results disagree with [12] in the squeezed limit, which used an interaction Hamiltonian 
where the leading order term is a total derivative, which complicates the calculations (see 
section 3). Upon comparing (4.27) with 

(C(r7,ki)B(r,,k2) • B(T7,k3)) = b'^'f'{2Trf6^''Hki + k2 + k3)Pc(^i)^i?(fc) (4.28) 
in the squeezed limit, we obtain the local type magnetic non-linearity parameter 

in agreement with the magnetic consistency relation [7]. 
5 Discussion and conclusions 

In a model where conformal invariance of a C/ (1) vector field is broken by a direct coupling of 
the vector field with the inflaton, we have calculated the three-point cross-correlation function 
of two powers of vector field with the primordial curvature perturbation, (C(ki)Aj(k2)Aj(k3)), 
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generated during inflation. We note in passing that the result can be generalized straight- 
forwardly to the case where C is replaced by an isocurvature field or a curvaton, under the 
identification in (3.16). 

In the case where the vector field is identified with the electromagnetic field, the strength 
of the cross-correlation is parametrized by the magnetic non-linearity parameter bj^L intro- 
duced in [7]. In the squeezed limit where the momentum of the curvature perturbation 
vanishes, the magnetic non-linearity parameter takes a local form 6^*^' = — (l/i/)(A//A/) , 
in agreement with the magnetic consistency relation [7]. On the other hand we find some 
disagreements with the results of [12], even in the squeezed limit. 

With the parametrization of the coupling of the form A((/)(r)) = Xi{t{t/ti)~'^" , the 
magnetic non-linearity parameter can be related to the spectral index of the magnetic field 
power spectrum ub = 4 — 2n for n > 0, giving the consistency relation 6^™^ = — 4. 
Thus in the most interesting case of a scale invariant magnetic field spectrum = 0, the 
non-linearity parameter is non- vanishing [7]. 

The full shape function is maximal in the flattened shape, where the momentum of the 
curvature perturbation is maximal and is enhanced by a logarithmic factor, in agreement 
with^ [8]. For length scales relevant for CMB, the logarithm will give an enhancement by 
a factor of 60. Thus, for a flat magnetic field power spectrum, the non-linearity parameter 
in the flattened limit becomes of order I^tvl] ~ 0(10^) depending on the scale. It will be 
interesting to understand if such a large b^i could be detectable. Previously it has been 
argued that the cross-correlation bispectrum of the metric perturbation with the magnetic 
fields can be observable through the combined survey of large scale structure and Faraday 
rotation measurements [8, 49]. We leave it for the future to explore this aspect in more 
details. 

On a different note, the identification of the vector field with the electromagnetic field is 
problematic from a theoretical point of view due to the strong coupling problem. Although 
we do not yet have a complete solution to the strong coupling problem, one can consider 
the model as a simple effective model for primordial inflationary magnetogenesis, assuming 
that a possible consistent UV complete model might not affect the generic phenomenological 
features [22]. 

Finally, we mention that the results obtained here might find applications beyond infla- 
tionary magnetogenesis. If the vector field is identified with a vector curvaton non-minimally 
coupled to the inflaton, then in a scenario where the CMB arises from a mixture of curvaton 
and inflation perturbations, the cross-correlation bispectrum of the curvature perturbation 
with the vector perturbations will be a new non-Gaussian signature in such models. In such 
a setup, there is no strong coupling problem, since the vector curvaton can be arbitrarily 
weakly coupled to the inflaton at the end of inflation. 
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Note added 



While writing up this paper, a related paper appeared [50] with a discussion of the effect 
on the CMB of the cross-correlation bispectra between the curvature perturbation and the 
vector modes. The authors of [50] use the same interaction Hamiltonian as in [12], where the 
leading order term is a total derivative (see our section 3 and the discussion). 

A Integrals 

Here, we evaluate the integrals In^ and In^ for different values of n. For n = 0, the integrals 
are given by 



For n = 1, the integrals are given by 
1 



^?^ = -7777317^(2^1+^2 + ^3) (A.2) 



■t 

{-kf - 2kl{k2 + k^) - 2ki{kl + k2k^ + kl) - {k2 + A;3)(fci + k2k3 + kj)) (A.3) 



For n = 3, the integrals are given by 
(1) _ 1 



Z 



{k2k^f/^k^t 

X [(A;i + k2f (3(-l + -f)kl - 3klk2 - k^) + {ki + fcs) ((-9 + 6j)kf - Qklk2 - 2kl) k^ 
+ (3(-3 + -f)kl - 6klk2 - 2kikl - 2kl) kl - 2{2kl + kik2 + kl)kl - 2{ki + ^2)^3 - A:| 
+ Skfk'flni-ktTi)] (A.4) 



j{2) ^ ^ 
3 2(^2^:3)7/2^2 

X [{45kfkf + (15(-1 + j)kl + 15(-3 + 2j)kf{k2 + hs) 

- 3kl{6k^ + 5(5 - 27)^2^3 + 6kl) - k\{k2 + A:3)((-49 + 307)(A;i + kl) + 10A;2A;3) 

+ Akl{k2 + hi){2,k'^ + 2k^k3 + Ak^kl + 2k2kl + 3/c|) 

+ kl{{A7 - 157)(/c| + k^) + (89 - 307)(/c^A;3 + ^2^1) + 6(16 - 57)A;|/c|) 

■~|~ (/^2 ~l~ ^2^3 ^2^3 ^2^3 ^2^3 ^2^3 ~i~ ^3) 

+ 2(^2 + fc3)(fef + A:f fea + A:|A;| + + fc^A:! + yfc2A:| + fcf)) 

+ lbklkj{kl -kl- kl) ln(-A:tr/)] (A.5) 
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